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Abstract 

We investigate the equation (u t + (f(u)) x ) x = f"(u)(u x ) 2 /2 where f(u) is 
a given smooth function. Typically f(u) = u 2 /2 or u 3 /3. This equation mod- 
els unidirectional and weakly nonlinear waves for the variational wave equation 
uu — c(u)(c(u)u x )x = which models some liquid crystals with a natural sinu- 
soidal c. The equation itself is also the Euler-Lagrange equation of a variational 
problem. Two natural classes of solutions can be associated with this equa- 
tion. A conservative solution will preserve its energy in time, while a dissipative 
weak solution loses energy at the time when singularities appear. Conservative 
solutions are globally defined, forward and backward in time, and preserve in- 
teresting geometric features, such as the Hamiltonian structure. On the other 
hand, dissipative solutions appear to be more natural from the physical point of 
view. 

We establish the well-posedness of the Cauchy problem within the class of 
conservative solutions, for initial data having finite energy and assuming that 
the flux function / has Lipschitz continuous second-order derivative. In the 
case where / is convex, the Cauchy problem is well-posed also within the class of 
dissipative solutions. However, when / is not convex, we show that the dissipative 
solutions do not depend continuously on the initial data. 
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1 Introduction 



A nonlinear variational wave equation whose wave speed is a sinusoidal function of the 
wave amplitude arises in the study of nematic liquid crystals. It is given by 

^-#9,(^)^ = 0, (1.1) 

with 

c 2 (» = asin 2 0) + (3cos 2 (ip), (1.2) 

where a and (3 are positive physical constants. We refer the reader to [11], [12], [14] 
for background information on the equation. Glassey, Hunter, and Zheng [10] have 
shown that singularities can form from smooth data for equation (1.1)-(1.2). Assuming 
that the wave speed c(-) is a monotone increasing function, the global existence of 
(dissipative) weak solutions has been established in [19], [20], [21], [23]. The general 
problem of the global existence and uniqueness of conservative solutions to the Cauchy 
problem of equation (1.1) will be addressed in a forthcoming paper [3]. 

The study of solutions to (1.1)- (1.2) consisting of a small-amplitude and high- 
frequency perturbation of a constant state has greatly contributed to the understanding 
of this equation [10], [19], [20], [21], [23]. Hunter and Saxton first studied these waves 
in [12]. Given a constant state a, these perturbed solutions take the form 

ip(t, x) = a + eu(et, x — c(a)t) + 0(e 2 ). 

Hunter and Saxton found that w(-, ■) satisfies 

(u t + u n u x ) x = l -nu n - x u 2 x (1.3) 

up to a scaling and reflection of the independent variables, assuming that a is such 
that S k \a) — 0, k — 1, 2, . . . n — 1, but c^ n \a) ^ 0, for an integer n > 1. In connection 
with our sinusoidal function c modeling nematic liquid crystals in (1.2), the relevant 
approximations in (1.3) are those with n — 1,2. The case n = 1 yields the first-order 
asymptotic equation 

(u t + uu x ) x = -u x , (1.4) 

for which existence and uniqueness of admissible conservative and dissipative weak 
solutions have both been established, see [13] and [16], [17], [18]. This equation is also 
an asymptotic equation of the Camassa-Holm equation [4], describing the motion of 
solitary waves in shallow water. For recent literature on the Camassa-Holm equation, 
we refer the reader to [5], [6], [7], [8], [9], [15] and in particular [2]. 
The case n = 2 yields the second-order asymptotic equation 

(u t + u 2 u x ) x = uu 2 x . (1.5) 

In [22] Zhang and Zheng established that dissipative solutions exist for (1.5) with BV 
data. In the analysis of (1.1), a major difficulty is concentration of energy at points 
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where d = 0, as in the example on p. 70 of [10]. We hope that investigation of 
singularities of the same type for the second-order asymptotic equation will be helpful 
toward the understanding of the original equation (1.1). 

Rather than (1.3), in the present paper we study a somewhat more general class of 
equations, having the form 

(«* + /(«)*)* = \f"(u)u 2 x . (1.6) 

Here u = u(t,x) is a scalar function defined for (t,x) G R+ x R + where M+ = [0, oof, 
and / is a C 2 function. More restrictions on / will be specified later. As initial and 
boundary data we take 

u(0,x)=u(x), u(t,0) = 0. (1.7) 

Integrating equation (1.6) w.r.t. x, we obtain 

u t + f\u) x = l -J*f\u)u 2 x dx. (1.8) 

It is now clear that, to make sense of this equation, we should require that the function 
u(t, •) be absolutely continuous with derivative u x (t, •) locally square integrable, for 
every fixed time t. Moreover, to satisfy the boundary condition at x — 0, one needs 
the nonnegativity of the characteristic speed at u — 0, namely 

/'(0) > 0. (1.9) 

The local integrability assumption u x (t, •) G L z 2 oc (R+) imposes a certain degree of 
regularity on the function u. Therefore, there is no need to consider weak solutions in 
distributional sense and a stronger concept of solution can be adopted. 

Definition 1.1 A function u = u(t,x)) defined on [0,T] x R + is a solution of the 

initial-boundary value problem (1.7) -(1-9) if the following holds. 

(i) The function u is locally Holder continuous w.r.t both variables t, x. The initial and 
boundary conditions (1.7) hold pointwise. For each time t, the map x \— > u(t,x) is 
absolutely continuous with u x (t, •) G L^ oc (IR + ). 

(ii) For any M > 0, consider the restriction of u to the interval x G [0,M]. Then the 
map 1 1— > u(t, •) G L 2 ([0, M]) is absolutely continuous and satisfies the equation 

ju{t, •) = -f(u)u x + \J* f'\u) u\ dx (1.10) 

for a.e. t G [0, T]. Here equality is understood in the sense of functions in L 2 ([0, M]). 

In spite of the regularity assumptions, the requirements contained in the above 
definition are still not enough to single out a unique solution. Let us consider a simple 
example. 
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Example 1. Consider the flux f(u) = u 2 and choose the initial data 



u(0, x) = 



-x, < x < 1, 
-1, X > 1. 



For t e [0, 1[ , a solution to (1.8) can be constructed by the method of characteristics, 
namely 



u(t, x) 



-x/(l-t), 0<x<(l-t) 2 , 



-(1-t), x>(l-tf. 

Notice that the norm of the gradient ||Ma;(^)||L°° blows up as t — > 1. For t — 1 we 
have = for all x > 0. At this stage, there are infinitely many ways to further 

prolong the solution. For example, we could set 

u(t,x) = t>lx>0. (1.11) 

Or else we could choose an arbitrary point b > 0, an arbitrary amount of energy k > 
and a time r > 1 and define 



u(t,x) =0 for 1 < t < t, 



while, for t > r, 



< x < 6, 

w(*,a;) = ^ (a; - b)/(t-r) b < x < k(t - t) 2 + b, 
■ (t-r) x>b + k(t-r) 2 . 

Among all these solutions, two in particular can be singled out. If we insist that the 
future configurations u(t,-) for t > 1 should be entirely determined only by the present 
configuration w(l, •), then the only reasonable choice is (1.11). On the other hand, if 
we look for solutions that satisfy the additional conservation equation 

(u 2 x ) t + (2uul) x = 0, 

the natural choice should be 

u(t x)-l x/{t ~^ 0< x <(t-l) 2 

To express the fact that at time t = 1 this solution is different from the null solution, 
in some way we should think its derivative u x as being not the zero function but the 
square root of a Dirac distribution at the origin. □ 

In the following, we say that a solution u = u(t,x) is dissipative if the family of 
absolutely continuous measures {/i( t ) ; t > 0} defined by dfi( t ) = u 2 (t) dx provides a 
measure-valued solution to 

w t + (f(u)w) x <0. (1.12) 
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More precisely, we require that 



<f)(t, •) dn {t ) 



< 



t2 



tl 



[<f> t (t,-) + <f> x (t,-)f'(u(t,-))]d m 



dt 



1.13) 



for every t 2 > h > and any function G C], > 0. 

On the other hand, to define a semigroup of conservative solutions we need to 
consider a domain V of couples (u, (i) where u : R + i— > R is an absolutely continuous 
function with square integrable derivative and /i is a nonnegative measure on M + . 
Decomposing fi = fi a + /i s as a sum of an absolutely continuous and a singular part 
(w.r.t. Lebesgue measure), we shall require that d\i a = u^dx. We say that a map 
t i— > (-u(t), /!(£)) G £> is a conservative solution of (1.7)-(1.9) if u is a solution 
according to Definition 1.1 and (1.13) is satisfied as an equality, for allt 2 > t x > and 

As mentioned earlier, Zhang and Zheng have established in [22] the finite-time 
singularity formation in smooth solutions and the global existence of a dissipative weak 
solution to (1.7)— (1.9) with initial data u(x) whose derivative is in BV, for f(u) = u 3 /3. 

In the present paper, we consider a general flux / with Lipschitz continuous second- 
order derivative such that f'(0) > 0. The initial data are chosen in the set of absolutely 
continuous functions m, with m(0) = and u x G L 2 (R + ). Our main results can be 
summarized as follows. 



1. A flow of conservative solutions can be globally defined, forward and backward in 
time (Theorem 3.1). The conservative solution of the initial-boundary value problem 
(1.7)— (1.9) is unique, provided that a suitable non-degeneracy condition is satisfied 
(Theorem 4.1). 

2. Assuming, in addition, that the flux / is convex, then there also exists a continuous 
semigroup of dissipative solutions. The dissipative solution of the initial-boundary 
value problem (1.7)— (1.9) is unique (see Theorem 5.1). 

3. If the flux / is not convex, the dissipative solutions do not depend continuously on 
the initial data, in general (see Example 2 in Section 6). 

Before proving the main results in Section 3, we briefly discuss the action principle 
and some admissibility conditions, whose aim is to identify a unique physically relevant 
solution to equations (1.7)— (1.9). 



2 Remarks on admissibility conditions 

The decay estimate 

u x (t,x) < 2/t 
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was used as an admissibility criterion for dissipative solutions of the first-order asymp- 
totic equation in [13], [16], [17] and [18]. We remark, however, that this is not appro- 
priate in connection with dissipative solutions of (1.6). Indeed, for a solution of the 
second-order asymptotic equation, the gradient u x can approach +oo as well as — oo. 

Another common criteria for selecting physically admissible solutions is by vanishing 
viscosity. One might conjecture that dissipative solutions are precisely the limits of 
vanishing viscosity approximations. We believe this is indeed the case when the flux 
function / is convex, see some proofs in [13] [15] for / = u 2 /2. On the other hand, 
when / is not convex, the dissipative solutions do not depend continuously on the 
initial data (see Section 6). We observe that the set of vanishing viscosity limits is 
closed, connected and depends on the initial data in an upper semi continuous way. 
Therefore, by a topological argument, the vanishing viscosity criterion cannot single 
out a unique limit, in general. 

Concerning the vanishing dispersion limit, numerical experiments performed with a 
convex / seem to indicate that vanishing dispersion selects the conservative solutions, 
see [13]. 

Next, we discuss the admissibility of solutions in terms of a variational principle. 
For all asymptotic equations (1.3) we have the action functionals 




In other words, the Euler-Lagrange equations satisfied by functions u that render sta- 
tionary the action A n are precisely the asymptotic equations (1.3). These can be 
derived from the nonlinear variational wave equation (1.1) 

i-cWW^J^O (2.2) 

by a perturbation argument. Notice that (2.2) is the Euler-Lagrange equation corre- 
sponding to the Lagrangean 

C = tf-c 2 ml (2.3) 
This arises often in physical models. For weakly nonlinear waves of the form 

xp = iJ + eu(r, 9) + e 2 v(r, 9) + 0(e 3 ) 

with 

t = et, 9 = x-c t, c = c(V'o) , 
assuming that c' = c'ijpo) 7^ we have 

iptt ~ c(ip)(c(ip)ip x ) x = -2c e 2 j(w r + c' o uu ) - ^c ^| + 0(e 3 ) . 

Moreover 

i) 2 t - c 2 {i))i) 2 x = -2c e 3 [u T u e + c' uu 2 ] + 0(e 4 ). 
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Therefore, u satisfies the first order asymptotic equation. The corresponding La- 
grangean, approximated to order 0(e 3 ), is —(u T ug + c' uu^). 

At first sight, one might hope that the physically relevant solutions to the equations 
(1.3) are those which maximize (or minimize) the action in (2.1). Unfortunately this is 
not the case, because the action A n is not coercive. For any smooth solution u of (1.3) 
one can find compactly supported perturbations u + ev which increase the value of A n , 
and others which decrease it. The extremality of the action thus cannot be used as a 
selective criterion. 



3 Conservative solutions 

We consider the evolution problem described by the equation 

1 r x 

u t + f( u ) x = - f"(u)u 2 x dx for all t > , x > 0, (3.1) 

2 Jo 

together with the boundary conditions 

u(t,0) = for all t > 0. (3.2) 

We assume that / G C 2 (IR) and 

/'(0)>0, |/»-/»|<L| U - V |, Viz, t) G R (3.3) 
for a constant L. 

One easily checks that every smooth solution satisfies the additional conservation 
law for the "energy" u 2 x1 namely 

(u 2 x ) t +[f(u)(ul)] x = 0. (3.4) 

It is therefore natural to seek a continuous flow associated with (3.1)-(3.2) which pre- 
serves the energy J °° u 2 (t , x) dx . However, Example 1 in the Introduction already 
pointed out a basic difficulty which one encounters while constructing a semigroup in 
the space Hl oc . Indeed, when the gradient u x blows up, all the energy is concentrated 
at a single point, so that the measure u x dx approaches a Dirac mass. 

Motivated by this example, to the equations (3.1)-(3.2) we will associate an evolu- 
tion semigroup on a domain T> defined as follows. An element of "D is a couple (u,fi), 
where u : M + i— > R is a continuous function with -u(O) = and whose distributional 
derivative u x lies in L 2 , while \i = fi a + fi s is a bounded nonnegative Radon measure 
on M + , whose absolutely continuous part (w.r.t. Lebesgue measure) satisfies 

dfx a = u 2 x dx. (3.5) 

In the following, on the family of Radon measures on IR + we consider the distance 



d(fi, ji) = sup 

<p 



ipdfx — / (p dfi 



(3.6) 
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where the supremum is taken over all smooth functions ip with \ip\ < 1, \ip x \ < 1. 

We recall that a semigroup S on a domain V is a map S : V x [0, oo[i— > D such 
that SoU> = u> and S s (S t w) = S s+t w for every s, t > and w <E V. 

Theorem 3.1 Assume that the flux function f satisfies condition (3.3). Then there 
exists a semigroup S : T> x [0, oo[i— > T> with the following properties. Calling t i— > 
S t (u,fi) = (u(t), ji(t)) the trajectory corresponding to an initial data (u,p) G £>, one 
has: 

(i) The function u = u(t,x) is locally Holder continuous in IR + x R + . provides a 
solution of (3.1)- (3. 2) in the sense of Definition 1.1 with initial condition 

u(0,x)=u(x). (3.7) 



(ii) The assignment 1 1— > /i( t ) provides a measure valued solution to the linear trans- 
port equation 

w t +\f'(u)w] x = 0, w(0) = /i. (3.8) 

Moreover, the singular part of the measure f"(u(t)) • fj,(t) vanishes at almost every time 
t > 0: 

f"(u(t))fi s {t) = 0, a. e. t. (3.9) 

(Hi) (Temporal continuity) For every M > 0, the above solution u and the corre- 
sponding measure ji satisfy the Lipschitz continuity property: 

i-M 

I \u(t, x) — u(s, x)\ dx < C\t — s\, (3.10) 
Jo 

where the constant C depends only on M , on the flux function f , and on the total 
energy /i(M+) < oo. 

(iv) (Continuous dependence on the initial data) Finally, consider a sequence of 
initial conditions (u n , Ji n ) G T> with u n — » u uniformly on bounded sets and d(ji n , fx) — > 
as n oo, for some (u, ft) G X>. Then the corresponding solutions satisfy 

u n (t,x) -» u(t,x) (3.11) 

uniformly for t, x in bounded sets, while 

<J(A*5),A*(t))->0 (3.12) 

for every t > 0. 
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Proof. We treat here the case where p, has compact support, say contained in the 
interval [0,i£], so that u is constant for x > R. The general case follows by an easy 
approximation argument. The proof will be given in several steps. 

1. Construction of the solution. Let an initial data (u, fi) G T> be given. Set 

£ = /i(IR+) < oo. On the semi-infinite strip {t > , £ G [0, £] } we construct a function 
U = U(t, £) by first setting 

U(0,O = U(O = u(m): (3-13) 

where 

y(0=mf{x>0; jl([0,x])>z} (3.14) 

for < £ < £, while 

y(0) = sup{x; /Z([0,x]) = 0}. (3.15) 

Observe that the map £ i— > y(£) is nondecreasing and left continuous, but it may well 
have upward jumps. The provision (3.15) makes it continuous at the point £ = 0. In 
any case, the composed function £ i— > u(y(^)) is always continuous. For positive times, 
the function {7 is then defined to be the solution of 

^0 = \l*f( U M)dri (3.16) 

with initial data (3.13). By the assumption of Lipschitz continuity of /", the function 
U can be obtained as the unique fixed point of a contractive transformation. Details 
will be given at the next step. 

Having constructed U(t,£), the characteristic curves are obtained by solving the 
equation 

y(o,0 = m, ^(t,0 = f'm,0). (s.ir) 

Explicitly: 

V(t,t) = m+ [ f'(U(r,0)dr. (3.18) 
Jo 

Notice that t i— > U(t, £) yields the values of our solution u along the characteristic 
curve 1 1— > y(t,£) starting at y(£)- A remarkable feature of equation (3.1) is that, if the 
energy is conserved, then these values can be determined in advance, before computing 
the actual position of the characteristic curve. The image of the mapping 

is now contained inside the graph of the desired solution u(t, •). More precisely, for any 
given (t, x) we define 

u(t,x)=U(t,£(t,x)), (3.19) 

where 

£(t,x) = sup{£; y(t,£) < x} . 
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Finally, at time t the corresponding measure is defined as the push-forward of 
the Lebesgue measure on [0, £] through the mapping £ i— > y(t,£). For each Borel set 
J C R + we thus define 

m (J) = meas G [0,£] ; j/(t,£) G j} . (3.20) 



2. A contractive transformation. Consider the space of continuous functions 
C([0, oo[ x[0, £]), with the equivalent weighted norm 



sup e- L ^\U(t,0\ , 
t>o, £e[o,£] 



(3.21) 



where L is a Lipschitz constant for the function /". The transformation U \— > TC/ is 
defined as ^ 

re/ (*, o = u(m) + lJ o J o f " {U{s ' v)) dv ds ' (3 ' 22) 

If now \\U — V\\* — 5, then 

f"(U(T,T,)) - r(V(T,r,)) < L\U(r, V ) - V(r,r,)\ < Le L ^5. 



For every t > and £ G [0, £ ] we thus have 



\(TU-TV)(t,0\<l " 



dr 



<-[ 

~ 2 J 



1 - 1 - 

s 2 



By the above inequality we conclude 



\\TU-TV\U<-\\U-V\U, 

proving the contractivity of the map T. By the contraction mapping theorem, it 
admits a unique fixed point U — U(t,£), defined on M + x [0, £ ] . In turn, the function 
u = u(t,x) and the measures fj,^ are well defined by (3.19)-(3.20). 

3. Absolute continuity. We prove here that the map £ i— > U(t,£) is absolutely 
continuous, for each t > 0. Indeed, consider first the case t — 0. Let with 
fc = 1, . . . , AT, be disjoint intervals contained in [0, £]. Assume that 

£l&-&l<e- 



Call Ii the set of indices k such that 

l^)-^(&) 



v(&) - !/(&) 
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and let 7 2 be the set of indices where the above quantity is > ^fe. Then 
keh keh 



while 

<^£&-&l<V£- 



fce/ 2 

Together, the two above inequalities yield 

AT 



fc=l 



proving the absolute continuity of the map £ i— > i7(0, £). 

For t > 0, the absolute continuity of U(t, •) follows from the absolute continuity of 
[7(0, •) together with (3.16). Indeed, 

if(f,o - < i^(o,o - u(o,o\ + sup irni . 

£ u 

As a consequence, the partial derivative U% exists at a.e. (t,£). By (3.16), it satisfies 
the evolution equation 

^-Udt,0 = \f"(U(t,0)- (3.23) 
On the other hand, the map £ i— > y(t,£) can be discontinuous. However, if 

lim y(i,f) = yi < y 2 = lim 2/(t, , 

then the function u(t, •) must be constant on the interval [2/1 , 2/2] - 

4. Measure transformations. To proceed, we first need to analyse the regular 
and the singular part of the push-forward of Lebesgue measure, under a continuous 
non-decreasing transformation. 

Lemma 1. Let U : [0, £] 1— > R fre absolutely continuous with square integrable deriva- 
tive. Let £ 1— > y(£) fre snc/i #ierf 

y(0 = I/(«)+ ful(()d(. (3.24) 

./0 
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For x G [y(0), y(C)] define the function u = u(x) implicitly by 

u(y(0) = ^(0 • (3-25) 
Let (j, be the push-forward of Lebesgue measure on [0, £] through the map y, i.e. 

fi(J) = meas{^e[0^}; y(£) E J} . (3.26) 

Then the absolutely continuous and the singular part of ji w.r.t. Lebesgue measure 
are respectively given by 

f i a (A) = meas{^e[0,^]; y(£)eA, U&) ^ 0} . (3.27) 

»'(A) = meas{te[0,£\; y(£)eA, U ( (t) = 0}. (3.28) 
In addition, on the set [y(0), y(Q] one has 

dfx a = u 2 x dx. (3.29) 

Viceversa, if both U and the map y are absolutely continuous and (3.25), (3.26), (3.29) 
are valid, then (3.24) must hold. 

Proof. By (3.24), the image of a set I C [0,£] 

7 £ = {££[(),£]; \U&)\<e} 

under the mapping £ i— > has Lebesgue measure 

meas(y(J)) = yVf(0^. 

It is thus clear that the singular part of /x is the push-forward of Lebesgue measure on 
the set 

7 = {££[(),£]; tf € (0 = O} 
Next, for any fixed e > take a measurable set J C [0, £] such that 

C/|(0>£ for all J. 

Then 

/ u 2 x (x) dx = [ lilt % ■ di = [ [Uz U7 2 ] 2 U\ ■ di = meas (J) . 
Jy{j) J j I dy] d£ J j 

Since e > was arbitrary, this proves (3.27), (3.29). To prove the last statement, 
assume (3.25), (3.26) and (3.29). Call 

J. = {ee[0,f]; y € (0>e}. 
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Observe that, for £ G J e . the chain rule yields 

«»(y(0)i/«(0 = ^(0- (3-30) 

For 0<a<5<£we now obtain 

/ u 2 x {x)dx= [ ^(y(0)^(0^ = meas([a,6]n J e ) (3.31) 

Since a, 6 were arbitrary, this implies 

%(0 = [^(y(O)] -1 (3-32) 

for £ G J £ . Together with (3.30) this yields 

«x(y(0) = ^ _1 (0, ife(0 = o?(0 (3-33) 

for all J £ . Since £ > is arbitrary, we conclude 

y(0=y(o)+ / € % (C)dC = y(o) + iim/ v € (C)dC = v(o) + /V|(C)rfC, 



proving (3.24). 



5. A class of regular solutions. Having constructed the trajectory t i— > (w(t, •), 
we still need to prove that it satisfies equation (3.1), coupled with the initial and 
boundary conditions (3.7), (3.2). We carry out the analysis first assuming that the 
map £ i — ?/(£) is absolutely continuous. At the end, this assumption will be removed. 
For each t > and £ G [0, £ ] define 

j/(t, 0) = j/(0) + t /' (0) , y(t, = y(t, 0) + f Ul(t, C) d( . (3.34) 

Jo 

By (3.23) this implies 

^Vdt, = J^f = f'(U(t, 0) tf € (f , • (3.35) 

We now check that the function y = y(t,£) defined at (3.34) coincides with the one 
defined at (3.18). Indeed, by the second part of Lemma 1, their derivatives y% coincide 
at time t — and satisfy the same equation (3.35). In particular, from (3.34) it is clear 
that the map t h- > y(t,£) is non-decreasing. In particular, characteristics never cross 
each other. 

We begin by observing that the boundary condition (3.2) is clearly satisfied, because 

u (t, 0) = U(t, 0) = 17(0, 0) + / U t {r,0)dT = 0. 

Jo 
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Moreover, the initial condition (3.7) holds because of the definitions (3.13)-(3.14). 

To check that the limit function u satisfies (3.1), fix a time t > 0. Since u(t,x) = 
for x G [0, 2/(t, 0)], in this region the equation (3.1) trivially holds. For almost every 
x G [?/(£, 0), y(t,£)], there exists a unique £ e [0, £] such that x = y(t,£). In this case, 
our construction yields 

i /•£ i rv{t,£) 

u t + f{u)u x = U t (t, = 3/ f"( U ^ 0) ^ = 2 X ' )} rf/i(t) • 

This implies (3.1), provided that we can show the identity of measures 

f"(u)u 2 x dx = f"(u)d m (3.36) 

for almost every time t > 0. We shall now work toward a proof of (3.36). 
Since the function u is continuous, by covering the open region 

{(M)£t + xR + ; f'WM))^} 

with countably many sets of the form 

T = [(t,x); te[h,t 2 }, xe[y(t,a), y(t,b)}} 

it suffices to prove the following statement. 
Assume that 

f"(u(t,x))>S>0 (t,x)eT. 

Then, for a.e. t G [ti, t 2 ], the restriction of the measure H( t ) to the interval [y(t, a), y(t, b)] 
is absolutely continuous w.r.t. Lebesgue measure and satisfies d/i^ = v? x dx. 
By construction, as long as U ranges in a region where /" > 5 we have 

|w.e>5- (3.37) 

Hence, for any e > 0, 

meas({(t,0 GT; |t^(i, £)| < e}) < y £ ■ 

Since e > here is arbitrary, we conclude that there exists a set of times Af of measure 
zero such that 

meas({£G [a, 6] ; C/ e (t,0 = 0}) =0 

for all times t ^ A/". By Lemma 1, t ^ jVthus implies that the restriction of fi( t ) to the 
interval [y(t, a), y(t, b)] is absolutely continuous w.r.t. Lebesgue measure. Furthermore, 
by (3.34), Lemma 1 shows that its density is dfj,^ = u 2 x (t) dx. This concludes the proof 
of (i) and (ii) in Theorem 3.1, at least in the case where the function £ 1— > is 
absolutely continuous. 
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6. General initial data. We now consider a more general initial data (u, p) G T>. 
In this case, the map £ h- > |/(£) is non-decreasing, left continuous but not necessarily 
continuous. Its distributional derivative is thus a measure, say D^y = a = a a + a s . By 
the assumptions, the absolutely continuous part satisfies 

da a = u\d£, 

so that 

Jo 

Consider a new initial condition (u*,£L*) defined by setting 

^(0 =y(0)+ f ul(()d(, u*(y*(0) = u(0 

Jo 

/T(J) =meas{£; y*(0 G J} . 

By construction, for this new initial data the function £ ^ y*(0,£) = is abso- 

lutely continuous. Hence, by the previous analysis, the corresponding function u*(t,x) 
provides a solution to the initial-boundary value problem (3.1)-(3.2) with initial data 
(u*,p,*). It is now easy to check that the function constructed in (3.16)-(3.19) for the 
original initial data u satisfies 

u(t, y(t,t) + o*([0,t\))=U(t,t). 

More precisely, 

u(t,x) = U(t,£) where £ = inf {C 5 y(t, C) + <t"([0, CD > x} . 

By the previous analysis, u* provides a solution. Hence the same is true of u. 

7. Continuity properties. Recall that £ = /x(M+) < oo is the total mass of each of 
the measures /X( t ). We have 

Tot.Var.{w(t, •) ; [0, M] } < sJ^M . 
Since u(t, 0) = 0, for any x G [0, M] we have 

\u(t,x)\ < yf^M. 
This implies the Lipschitz continuity property w.r.t. time: 

\u(t,x) - u(s,x)\dx <\t-s\- {sup • y^M + ^ • sup (3.38) 

where both sup are taken over \u> 
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Next, consider a convergent sequence of initial data (u m , fi m ) m > . The assumption 
of Theorem 3.1 implies that the corresponding functions U n satisfy 

tno - mo 

uniformly on [0, £]. Therefore U m (t, £) — > U (t, £) uniformly on the domain [0, T] x [0, £] , 
for any T > 0. In turn, this implies the convergence (3.11)-(3.12). 

8. Holder continuity. We show that u(t,x) is Holder continuous locally in (t,x). 
First we know by Sobolev embedding that u is Holder continuous in x for each fixed 
time t with exponent a = 1/2. In the time direction, we know that the derivative 
of u along a characteristic is bounded, thus u is Lipschitz continuous in time along a 
characteristic. The characteristic speed is u which is locally bounded, thus the distance 
traveled in the x direction is order one of time. Combining the two parts, we conclude 
that u is Holder continuous locally in both space and time. 
This completes the proof of Theorem 3.1. 

Remark. The previous construction of solutions to (3.1)-(3.2) works equally well for 
negative times. The semigroup S can thus be extended to a group ^ : V x ffi. i— > V. 



4 Characterization of semigroup trajectories 

In the previous section, a solution u to the initial-boundary value problem (3.1)-(3.2), 
(3.7), was obtained as the fixed point of a contractive transformation. Hence, any other 
solution which provides a fixed point to the same transformation necessarily coincides 
with u. A straightforward uniqueness result can be stated as follows. 

Theorem 4.1 Assume that f satisfies (3.3). Consider a function u = u(t,x) and a 
family of measures /i^ satisfying (i) and (ii) in Theorem 3.1. Moreover, calling 

y(t,0 = inf{x>0 ; m ([0, x}) > ^} , (4.1) 

U(t,0=u(t,y(t,0), (4.2) 

assume that for a. e. £ the map t i— > U(t, £) is absolutely continuous and satisfies the 
differential equation (3.16). Then one has the identity 

(u(t),fi(t)) = S t (u,ji). (4.3) 
In particular, the solution which satisfies the above conditions is unique. 

We conjecture that a uniqueness result remains valid even without the assumption 
(3.16) on the corresponding function U. The basic ingredient toward a uniqueness 
result is the assumption 

f"(u)d^(t) = f"(u)ul(t)dx. (4.4) 



16 



for a.e. t. We now show that this is indeed the case under the additional condition 
f" > o. 

Theorem 4.2 In addition to assumption (3.3), let /"(•) > 0. Consider a function 
u = u(t,x) and a family of measures /i^ satisfying (i) and (ii) in Theorem 3.1. Then 
identity (4-3) holds. 

Indeed, observe that the flow on I/QO, £]) generated by the evolution equation 
(3.16) is Lipschitz continuous w.r.t. time and to the initial data. Adopting a semigroup 
notation, call t i— > V(t) — S t V the trajectory corresponding to the initial data V G 
L 1 ([0, £]). Since the couple (u(t), fiu\) can be entirely recovered from the function 
U(t, •) and the initial mapping £ i— > to prove uniqueness, it thus suffices to show 
that 



lim — 

h^0+ h 



U(t + h,0-(ShU(t))(0 d£ = (4.5) 



for almost every time t > (see Theorem 2.9 in [1]). Since /" > 0, our assumption 
implies that the singular part of /i(t) vanishes at a.e. t. Choose a time t where /z^ = 0. 
Then 

U^O^O fora.e. £ G [0,£] . (4.6) 
Consider the map £ h- > y(t,£). Since 

2 J AC tt d y 1 

u x dy = d£, Ut = u x - — = — , 

by (4.6) the pre-image of a set of measure zero through the map £ i— > £) has measure 
zero. 

If now if = u(t,x) is differentiable at the point (t, £)), we have the identity 

ft/(U)= [«t + /'(nK](U(U)) (4?) 

= \ /» ^ *) <** = I Jo m*> ^)) ^ • 

Observing that u(t, ■) is differentiable at a.e. x, we conclude that (4.7) holds at 
a.e. £ E [0, £ ]. In turn, this implies (4.5), proving the theorem. 



□ 

Notice how the condition on the vanishing of the singular part is essential to ensure 
uniqueness. Otherwise, in Example 1 the solution u(t,x) = for t > 1, with /i( t ) 
containing a unit mass at the origin, would satisfy all the other requirements of the 
theorem. 
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5 A semigroup of dissipative solutions 



Next, we examine dissipative solutions. A major difference with the conservative case 
is that here the Cauchy problem is well-posed if the flux function / is strictly convex, 
but ill posed otherwise, as shown in the next section. 

In this section, our main concern will be the construction of a semigroup of dissi- 
pative solutions under the additional assumption that /" > 0. As domain V of our 
semigroup we choose the space 

V = : 1R + i— > K. , u is absolutely continuous, u(0) = , u x G L 2 | . 

Theorem 5.1 Assume that the flux function f satisfies (3.3) and f" > 0. Then 
there exists a semigroup S : V x [0, oo[i— > V with the following properties. Calling 
t i— > u(t) = S t u the trajectory corresponding to an initial data u G V, one has: 

(i) The function u = u(t,x) is Holder continuous. It provides a solution of (3.1)- 
(3.2) with initial condition u(0, x) = u(x). 

(ii) For every M > 0, the above solution u satisfies the Lipschitz continuity property 
in time: 

I \u(t, x) — u(s, x)\ dx < C\t — s|, (5.1) 
Jo 



(Hi) Given a sequence of initial conditions u n G V, assume that 

K - ^x||l 2 ([o,m]) — »• 
for every M > 0. Then the corresponding solutions satisfy 

u n (t,x) -> u(t,x) (5.2) 

uniformly for t, x in bounded sets. 

Proof. Consider an initial condition u G V. For simplicity, we again assume that u 
is constant outside a bounded interval, say [0,R]. The general case follows from an 
approximation argument. 

To construct the corresponding trajectory we begin by setting 

f R 

f = / \u 2 x (x)\dx. 
Jo 

Then we define the initial data 

U{£) = am)) , 
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where 

y(0 = inf {x > ; J u 2 x (x) dx > ^ . (5.3) 

By the analysis in Section 3, the map £ i— > £/(£) is absolutely continuous, hence its 
derivative 

m = 

is a well defined function in L 1 ([0,£]). 
Define the subset 

J" = {£G [(),£]; Z(0<o}. 

Let L be a Lipschitz constant for /". On the space of continuous functions F : R + h- > 
L x ([0, £]) with weighted norm 

11^11* = sup e-^||y(0IUi, 

t 

we now define a continuous transformation F i— > TF as follows. 

TY(t, = Z(0 + jT Hv, Y{s, 77)) ^ ^ , (5-4) 

where 

$(77, F) = min{F,0} if r] E J~ , 

$(t7,f) = f if v e[0,t]\J-. 

To check that T is a strict contraction, assume that ||F — F||* = k, so that 



/ |F(t,0-F(t,OI^<^e L f 
./o 



for all t > 0. Then for every s > 



f $(77, F(s, 77)) — $(77, F(s, 77)) ^< K e^, 
./o 



and therefore 

. _~ /* /*Lk 



J \(TY-TY)(t,t)\dt<J J ^e L ^dsdi< K -e^. 



By the definition of our weighted norm, this implies 

\\TY — TFIL < - IIF-FI 

11 11* - 2 11 l 
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Let now Y = Y(t,£) be the unique fixed point of the transformation T. Then one 
easily checks that the function 

Z(t,£) = Y(t,0 if ZtJ-, 

Z(t,0 = min{F(t,0, 0} if ie J', 
provides a solution to the equations 

z(o, = ^(0, 



dz_ 

~dt 



(t, = \f" (J Z(t, 77) dr^ if + for all s G [0, t] , 
— (f , = if Z(s, = for some s G [0, t] . 



In turn, we can now define 

./o 

and the characteristic curves 



y(i,0 = y(0+ f f'(u(s,t))ds. 

Jo 



In a similar way as in Section 3, the dissipative solution u can now be obtained by 
setting 

u(t,x) = U(t,Z{t,x)), 

where 

£(t,x) = sup{£; y(t,g) < x} . 
To see why this construction actually yields a solution to (3.1), consider first the case 

where the map £ i— > is absolutely continuous. Then y^(0,£) = [/^(O, £) = Z 2 (0,£)- 
Since 

= /» ^ = f"(u) z = ^ 

for all t, £ we deduce the identity 

ydt,0 = Z%0 = U*(t,0. 
Moreover, (3.34) again holds. As in the proof of Theorem 3.1, we obtain the relations 

For almost every x G 0) , y(t, £)] , if x = y(t,£), then 

[u t + f'(u)u x ](t,x) = -u(t,y(t,£)) = -U(t,Z) = J -^Z(t,ri)dri 
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= \J rmt, v)) di = \J* /"(«(*. v)) <^ y) d y ■ 

The second identity in (5.5) was used here to change the variable of integration. 

The extension to the case of general initial data, where the map £ i— > y(t,£) is not 
necessarily absolutely continuous, is carried out as in the earlier proof of Theorem 3.1. 
We skip the details. 



6 Instability of dissipative solutions for non-convex 
flux 

In this section, we show that if the convexity assumption /" > is dropped, then the 
Cauchy problem for the equation (3.1)-(3.2) is ill posed, in general. 

Example 2. Consider the flux function f(u) = u 3 . Let U = U(t,£) be a solution of 
(3.16), with £ G [0, 3], such that at some time t > there holds 

re, e e [o, i] , 

U(t ,0 = { 2-£, £g [1,2], 
[0, £e[2,3]. 

Consider first the conservative solution u = u(t,x). This is well defined forward and 
backward in time. At time t = to, its explicit values are 

(x, are [0,1], 

u(to,x) = <2-x, x G [1,2], 

[o, x > 2 

while a unit mass is concentrated at the point x — 2. Assuming £ sufficiently small, 
we have 

= / 3C/(i, 77) rf?7 > 
for all t G [0, i ] and £ G ]0, 3]. Hence 

J^(*,0=3C/<0, f/ 5 (t,O<0 for tG[0,t [, 2<£<3. 

Next, consider a dissipative solution t> coinciding with u at time i = 0. This means 
v(0, x) = u(0, x) = U(0, for x = y(0, . (6.1) 

We recall that 

= / Ul(t,r))dr). 
Jo 
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Clearly, v will still coincide with u as long as its gradient remains bounded (equivalently, 
as long as remains bounded away from zero). On the other hand, for t > t , the 
dissipative solution v = v(t,x) coincides with the conservative one only on the interval 
where x < y(t, 2), while v is constant for x > y(t, 2). In other words, 

v(t, x) — u(t, x) if te[0,£o]i 

l± \ (u(t,x) 0<x<y(t,2), ± u ^ , 

v(t,x) = < )' L ~ /7n\ if t G t , 2t . 

v ; \u(t,y(t,2)) x>y(t,2), 0J 
Energy dissipation occurs at time t — to, namely 

3 te[0,t 



DO 



v x (t, x) dx 



2 f > to- 



Next, consider a family of perturbed initial conditions, say 

U £ (0,0 = U(0,0+eHO, 

where is a non-negative smooth function, whose support is contained in [0, 1]. Since 
U I— > f"(U) = 6U is a monotone increasing function, by a comparison argument from 
(3.16) we deduce 

u%t,o>u(t,o 

for all e, t > 0, £ G [0, 3]. In fact, for a nontrivial we can assume a strict inequality: 

C/ £ (t,0>C/(t,0 *>0, £e[2,3]. 

For 2 < £ < 3 we now use the relations 

g-ui(t,0 = 3c/ £ (t,0 > 3c/(t,o = ^(U) , ^|(o,0 = ^(o,0 , 

and deduce 

UI(t,O>Uz(t,O>0 te[o,t ]. 

Moreover, for t > t and 2 < £ < 3 one has 

^C/|(t,O = 3C/ £ (t,O>3C/(t,O>0. 

Therefore, for each e > 0, the quantity £/|(i,£) is still strictly positive at time t = t 
and increases afterwards. It thus remains uniformly bounded away from zero. 

Since u x = U^ 1 , the above implies that, for any fixed e > 0, the corresponding 
conservative solution u e = u £ (t,x) has a uniformly bounded gradient. The dissipative 
solution thus coincides with the conservative one. As e — > 0, at time t — our 
construction yields 

IK(o) - «(o) || c o -> o , ||<(o) - « x (o) || L2 -> o . 
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However, when t > t and x > y(t, 2) the previous analysis yields 

lim u 6 (t, x) = u(t, x) 7^ v{t, x) , 

e^0+ 

where u,v are respectively the conservative and the dissipative solutions of (3.1)-(3.2), 
with the same initial data (6.1). The example proves that dissipative solutions do not 
depend continuously on the inital data. 

Remark. The previous example also shows that the family of dissipative solutions 
may not be closed. Since the set of solutions which are limits of vanishing viscosity 
approximations is closed and connected, we see that this set cannot coincide with the 
set of dissipative solutions. 
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